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O ■ Abstract 
O 

Q^. The displacement noise in the test mass mirrors of interferometric gravitational wave detectors 

5-H ' 

is proportional to their elastic dissipation at the observation frequencies. In this paper, we analyze 



one fundamental source of dissipation in thin coatings, thermoelastic damping associated with the 
dissimilar thermal and elastic properties of the film and the substrate. We obtain expressions for 
the thermoelastic dissipation factor necessary to interpret resonant loss measurements, and for the 
spectral density of displacement noise imposed on a Gaussian beam reflected from the face of a 
coated mass. The predicted size of these effects is large enough to affect the interpretation of loss 
measurements, and to influence design choices in advanced gravitational wave detectors. 
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I. INTRODUCTION 



While recent results indicate that the elastic losses available in bulk materials such as 
silica and sapphire are adequate to achieve the design goals for displacement noise in next 
generation interferometric gravitational-wave detectors, the losses associated with the mul- 
tilayer dielectric mirrors deposited on the faces of the mass are large enough to contribute 
significantly to the total noise of the system. The origin of these coating losses is not yet 
clear. In this paper, we investigate an intrinsic dissipation mechanism, thermoelastic effects 
associated with a thin film on a bulk substrate. Thermoelastic dissipation has been known 
since Zener's work in the 1930's. In homogeneous solids, it is associated with the irre- 
versible flow of heat driven by temperature gradients associated with strain gradients in the 
solid. These effects lead to the well known result for the damping of fiexural vibrations in 
a thin beam, where the heat flows from the side of the beam in compression to the side 
in tension. P] The maximum dissipation (/'max in this case is a function only of the material 
properties and not the beam dimensions, 0max = ~ Ea'^T/C, where E is the Young's 
modulus, a is the thermal expansion coefficient, C is the volumetric heat capacity, and T is 
the background temperature. The dissipation peak occurs at a frequency c<jp ~ 1/r where 
r ~ PC/k is the thermal diffusion time through a beam of thickness / given the thermal 
conductivity k. 

More recently, thermoelastic dissipation in homogeneous test masses has been analyzed 
as a source of noise in gravitational wave detectors, 01 where the heat can be viewed as 
diffusing radially in the compressed region associated with a Gaussian-distributed pressure 
field on the surface of the mass, as is used in Levin's method for analyzing displacement 
noise. P] Characteristic of both these examples is the presence of a nonuniform strain field, 
necessary to create a temperature gradient to drive the thermal diffusion in a homogeneous 
medium. 

In an inhomogeneous body, temperature gradients can be generated in a uniform strain 
field, so that thermoelastic dissipation can be expected even in the absence of stress or strain 
gradients. For the case of a coated test mass, if the thermoelastic properties of the film are 
different from those in the substrate, we can anticipate that thermal diffusion and hence 
thermoelastic dissipation will occur. An estimate of the size of the effect can be obtained 
by comparison with the fiexural damping of a thin beam. Replacing the thermal expansion 
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coefficient by the difference between these coefficients in the ffim and substrate, and assuming 
for simphcity that the other pertinent material parameters are the same, we have 0max ~ 
E{af — asYT jC , and again expect the peak response to occur for frequencies corresponding 
to the thermal diffusion time through the film. If we consider a film with the thermal 
expansion coefficient of alumina on a substrate with the thermal expansion of silica, and take 
the other parameters to be those of silica we find 0max ~ 3 x 10"'^, comparable to the elastic 
losses measured in optical coatings. 0, 0] For a 5-micron-thick film with these properties the 
dissipation maximum occurs at ~ 5 kHz, corresponding to the thermal diffusion time through 
the film. This frequency is in the range typically sampled by mechanical loss measurements, 
and not far from the frequency band of interest for gravitational wave detection. 

It appears that a more quantitative investigation of these effects is necessary to evaluate 
their implications for characterization of test masses as well as for gravitational wave detec- 
tors themselves. We consider in this paper two questions associated with the thermoelastic 
mechanism: what is the effective dissipation in the situation characteristic of resonator mea- 
surements of elastic loss, and what is the expected spectral density of displacement noise in 
the situation characteristic of a test mass in a gravitational wave detector. 

The key results from the analyses contained here are: (a) the derivation of an expres- 
sion for 0tot, II , the thermoelastic dissipation expected in a coated test mass undergoing 
deformations of the type expected in mechanical loss measurements and (b) the derivation 
of an expression for Sx{f), the power spectral density of thermoelastic displacement noise 
associated with the dielectric mirror coating on a test-mass substrate. 

A recent independent calculation of the spectral density of thermal noise in the low fre- 
quency limit agrees, for the case where the thermoelastic properties of the film and substrate 
other than the thermal expansion coefficient are identical, with the results given here taken 
to that same limit. 0, ^• 

In reference the difference between the expressions for thermoelastic thermal noise 
presented here and in for cases where the elastic properties of the film and substrate 
differ is noted. After discussions with the authors this difference has been resolved in favor 
of the results presented here. 
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II. SKETCH OF THE CALCULATION 



In this section, we sketch the calculation of the thermoelastic dissipation and the displace- 
ment noise in an inhomogeneous medium. Details of the calculation are given in section HVI 
and appendices. 

There are three steps to calculating the thermoelastic loss in the coating: 

1. Obtain the oscillatory thermal field associated with the zeroth-order elastic fields, 

2. Calculate the complex first-order elastic fields generated by the spatially varying os- 
cillatory thermal field, and then 

3. Calculate the power dissipated by the interaction of the zeroth- and first-order elastic 
fields. 

Throughout we will consider only linear thermoelasticity, retaining terms up to first order 
in the oscillatory thermal field. Therefore the stored energy can be taken to be proportional 
to the square of zeroth-order elastic fields, while the imaginary part of the product of the 
zeroth-order elastic fields and the elastic fields induced by the thermal wave represent the 
relevant average dissipated power. Zeroth and first order quantities are denoted by subscripts 
0, 1, respectively. 

The geometry we consider consists of a film of thickness / and thermal expansion coef- 
ficient a/ on a substrate with as whose thermophysical properties are all possibly different 
from those of the film. We take the surface normal to be in the —z direction, and the surface 
located at z = 0, so that the film extends from 2; = to 2; = /, and the substrate from z = I 
to z = 00. 

To simplify the analysis, we assume that the multilayer film can be approximated as a 
uniform film with appropriately averaged properties, and assume that the film is thin enough 
and the thermal diffusion length at the frequencies of interest short enough compared to any 
relevant transverse dimension (e.g. the dimensions of the object itself, or the radius of 
the Gaussian beam interrogating its surface) that only the thermal diffusion normal to the 
surface of the mass need be considered. Since the thermal diffusion lengths for frequencies 
of interest are on the order of or longer than the total film thickness, the description of the 
film in terms of its average properties appears reasonable, but we also consider the case of 
a film whose thermal expansion (but no other parameters) varies periodically through its 
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thickness, as a simple model to explore any unexpected effects that might arise from the 
neglected microstructure of the film. 

The point of departure for the calculation is the thermal diffusion equation, driven by a 
thermoelastic source term, which for the assumed one-dimensional heat flow can be cast in 
the form 3] 

where Oj{z) is the time- varying temperature with exp{iujt) time dependence assumed, kj = 
kj/Cj is the thermal diffusivity, Ej is the Young's modulus, uj is the Poisson ratio, T is the 
background temperature, Cj is the heat capacity per unit volume, eo,Mj is the zeroth-order 
i-polarized compressional strain, and j = f,s indicates quantities evaluated in the film and 
the substrate, respectively. 

To obtain the source term, we need the zeroth-order compressional strains. Different 
combinations of zeroth-order elastic fields are relevant in different situations. We will as- 
sume that the transverse variation of the zeroth-order elastic fields is slow compared to the 
thickness of the film, so their only possible variation is in the z-direction, and that that 
variation results only from the possible discontinuity of the elastic properties at the film- 
substrate boundary. Note that this statement regarding the ^-dependence applies only to 
the zeroth-order elastic fields; as we will see, the thermal fields and the first-order elastic 
fields they generate have a z-dependence that arises from the propagation of the oscillatory 
thermal wave itself. Under these assumptions, we can specify the zeroth-order fields in terms 
of three quantities that, due to the elastic boundary conditions, do not vary over the length 
scales relevant to this problem: the axial stress o"o = cro,zz, the symmetric combination of 
in-plane strains (the dilation) Eq = {eq^xx + eQ^yy)/2, and the antisymmetric combination of 
in-plane strains £o,xx ~ ^o,yy of the other components of the zeroth-order elastic fields 
can be derived from these three, as is established in Appendix ^ The anti-symmetric com- 
bination of strains, which is a pure shear along axes rotated 7r/4 to x and y, does not interact 
thermoelastically, and can be neglected in the remainder of the analysis, as is established in 
more detail in section llV(^ 11 

Given these zeroth-order elastic fields, we can evaluate the source term in Eq. ^and solve 
for the oscillatory thermal wave, 0{z), as discussed in section ITV Al and Appendix |B| This 
thermal wave, in turn, generates a first-order elastic field, with compressional components 
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Ei^ii and (Ti^ii. The thermoelastic coupling enters into the formulation through the elastic 
equilibrium equations and modified Hooke's law, which can be cast for the one-dimensional 
case considered here from Eqs. 7.8 and 6.2 of [isf, 

^ [ei,xx,j + ei,yy,j + 2(1 - - 2(1 + Uj)aej] = (2) 



and 



1 + 



EjUjOj 



(3) 



where j = f,s represents fields and material properties in the film and substrate, respectively. 
With the thermal fields obtained by solving Eq. Q we can obtain from Eqs. |21 and |21 the 
first-order elastic fields, as given in Appendix O 

The rate at which work is done per unit volume in a deformed body is in general given 

by 

P = ""^'^ 

and the dissipated power density by 

Pdiss ^ -- - ^ Im [aliiSi^ii + alii€o^ii\ . (4) 
1=1 

where the second form is specified to our problem with sinusoidal fields and only longitudinal 
strains, and takes into account that the zeroth-order fields are real, so that dissipation first 
appears in the product of first- and zeroth-order fields. Integrating the dissipated power 
density over < 2; < oo, we obtain the dissipated power per unit area, given in section 
IIVBI While this is the essential quantity of interest, it is convenient for comparison to 
experimental measurements of elastic Q to define an effective dissipation factor, (p. Since 
the thermoelastic dissipation is nonlocal, the choice of stored energy with which to make 
such a definition is somewhat arbitrary. A reasonable choice, and the one that results in 
a value for directly comparable to that derived from experimental results, is the elastic 
energy stored in the film. With this choice, as described in Section HV C 1\ we find that for an 
elastic field with specified in-plane strain and vanishing axial stress, as would be appropriate 
for a measurement of the elastic Q of a mass coated with a uniform film of thickness /, the 
total thermoelastic loss, (j)tot (Eq- l55|) is well approximated by 



'^_asE,{l-uj)Cf'' 



afEf{l-v,)Cs\ 



where the frequency dependence is contained in the function g^u) defined by 

, , , r 1 sinh(y^) ] 

g{u;) = lm --=_—-==-————-_ , (6) 

JtujTf cosh[JtujTf) + nsmh[Jtujrf) 

where u = 271 f,Tf = I'^/nf = t^Cf/kj is the thermal diffusion time across the film, and 



R = ^kfCf/ksCs, with kj and Cj the thermal conductivity and volumetric heat capacity, 
respectively. Eqs. El and IHl (or the form for a multilayer in Eq. ^ constitute one of the 
two key results of this paper. The frequency dependence represented by g{uj) is discussed at 
length in section IIV C 31 The quantitative implications for measurements of thermoelastic 
dissipation in several material systems are discussed in section UllAI 

These numerical results indicate that thermoelastic losses associated with the coating are 
comparable to those obtained in experimental measurements of elastic loss, which suggests 
that their contribution to the total displacement noise budget for a test mass could be 
significant. While one could form an expression for the thermal noise imposed on a Gaussian 
beam interrogating a coated test mass by inserting (pi^ y from Eq. and a corresponding one 
for 4>i,± from Eq. I^into one of the expressions developed for thermal noise in coated test 
masses 0, 0] , the result would only be approximate because of the thermoelastic coupling 
between in-plane and normal strains. A direct calculation of the thermal noise can instead be 
carried out using Levin's formulation, calculating the power dissipated by the thermoelastic 
mechanism when a pressure field with the same radial distribution as the optical intensity 
field is applied to the coated mass. This calculation, using the zeroth-order elastic fields 
obtained in Isl , is carried out for arbitrary frequency in section IIVDI We find 

S.if)df = ^-^^{l + u^rA'gHdf 



TV^f W'^ C 

where is a dimensionless positive-definite combination of material constants that vanishes 
when the film and substrate are identical. 



1 + Uf , , Er' 



2 



-1 , (8) 



g{uj) is the same frequency dependence as defined in Eq. [HI and the second form holds for low 
frequencies obeying u < 1/r/. Note that the results in section UlI Bl show that the limiting 
form must be used in the gravitational-wave-detection band only with caution. Eq. [7| (or 



the corresponding Eq. QHlfor a multilayer) is the other key result of this paper; quantitative 
implications for several plausible mass/coating combinations are presented in Section Fill Bl 
While the results for 0/^ y and S^if) in Eqs.|niand[7|are calculated for a uniform film, most 
optical coatings of course will consist of a large number of layers. In such a multilayered 
coating, there are two thermoelastic dissipation peaks, one at a frequency related to the 
thermal diffusion time through the entire film, and one at a frequency related to the thermal 
diffusion time through an individual layer. These interlayer effects are investigated in section 
HVT^ It is seen there together with section IIV C 31 that for problems of interest, the 
thermoelastic effects are dominated by contributions from the thermal diffusion through the 
film, so that a description of the multilayer in terms of a set of averaged properties appears 
appropriate. The subtleties of the averaging process are investigated in Appendix El where 
it is seen that the average of various products of material quantities is required in addition 
to the average of the quantities themselves. Specializing to a periodic multilayer of total 
thickness / with alternating layers of materials a and b in thicknesses of da and d^, the 
result for 0; j| from Eq. |321 

2 

g{u;) (9) 



2CfT 



1 / Ea \ 1 E,a, 



Cp \1 — Z//avg C*s 1 



where the frequency dependence is contained in the same function g{uj) as defined in Eq. 
with Tf — > Tp where Tp = I'^/np, and the result for S^if) from Eq. [TH 



^ + ^sY^^=l^\ (10) 



and Eq. ESI 



C, ( a 



l + v , , E 

l + Us Es 



2 



^ 2asCp V 1 - 

can be stated in terms of an averaging operator defined in Eq. ID II as 



-1 , (11) 

avg J 



(-^)avg = -1 + ^—f^b , (12) 



and volume-averaged material properties Cp and Kp defined in Eq. ID 51 The second form 
of Eq. El is, like the second form of Eq. [3 a low frequency limit valid foru; < 1/t/. 
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III. NUMERICAL RESULTS 



A. Numerical results for dissipation in Q measurements 

The mechanical loss factors (or equivalently Q factors) of dielectric coatings applied to test 
mass substrates may be obtained experimentally. In a typical measurement of this type, a 
subset of the vibrational resonant modes of a coated substrate are individually excited above 
some background level and the decay of the amplitude of vibration of the face of the sample 
measured as a function of time. From this measurement, and a model of the distribution 
of the stored energy in each mode of a coated sample, the mechanical loss factors of the 
dielectric coating at each of the resonant frequencies of the sample may be obtained. 

As part of the loss measurement process described above, the coated samples experience 
periodic strains. If there exists a difference in the thermoelastic properties of the dielectric 
coating and the substrate, then as shown in this paper, there will be thermoelastic dissi- 
pation. Eq. ini can be used to calculate the thermoelastic dissipation in a coating both at 
the frequencies typical of mechanical loss measurements, and at frequencies of interest for 
gravitational- wave detection. 

The expected thermoelastic loss associated with a coating on a substrate is a direct func- 
tion of the material parameters for the particular substrates and coatings chosen. Current 
interferometric detectors use fused-silica substrates with coatings formed from alternating 
layers of Si02 (refractive index n = 1.45) and Ta205 (refractive index n = 2.03), each layer 
being of A/4 optical thickness, with A = 1064nm. The mirrors in future upgrades to current 
detectors are expected to have substrates of either fused silica or sapphire. The choice of 
appropriate m.ror coating .natena. is a subject of ongoing study Q. w.tK the two 
coatings currently under most intense study being alternating layers of Si02 and Ta205, or 
AI2O3 (refractive index n = 1.63) and Ta205. 

To estimate the expected level of thermo-elastic loss for the mirror/substrate coating 
combinations above, Eq. IHl was used. In each case a coating thickness equivalent to thirty 
alternating quarter-lambda layers of the coating materials was chosen. The numerical values 
used for the properties of the mirror substrates are shown in the tables below. 

Choosing appropriate material parameters for the multi-layer ion-beam-sputtered dielec- 
tric coatings is made more difficult by the fact that thermo-physical properties of these 
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TABLE I: Fused silica substrate properties 




as = 5.1 X 10"^ K"^ 

Es = 7.2 X 10^° 

Cs = 746 Jkg-^K-i X 2200 kgm-^ 

= 1.64 X 10^ 
ks = 1.38 Wm~i 
i^s = 0.17 



TABLE IL Sapphire substrate properties 




as 


= 5.4 X 10-6 


Es 


= 4 X 10" Nm-2 


Cs 


= 777 Jkg-^K-i X 3980 kgm 




= 3.09 X 10*^ 


ks 


= 33 Wm-i 



Us = 0.23 



types of coatings are not well characterized. Absent better information, the properties of 
the amorphous Si02 and AI2O3 present in the films were assumed to be the same as the 
bulk values for amorphous fused silica and crystalline sapphire. The numerical values used 
for the properties of Ta205 are summarized in table ITTTl 

TABLE IIL Properties used for Ta205 in thin film form. 



a 


= 3.6 X 10"*^ 


E 


= 1.4 X 10"Nm-2 


C 


= 306 Jkg-iR-i X 6850 kgm 




= 2.1 X 10^ JK-im-3 


k 


= 33 Wm-iR-i 


V 


= 0.23 



Values for Young's modulus and density of Ta205 in thin film form were taken from 
reference The coefficient of thermal expansion for Ta205 film was taken from and 
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FIG. 1: Calculated frequency dependence of the thermoelastic losses of fused silica substrates with 
dielectric multilayer coatings formed from thirty alternating quarter-lambda layers of (a) Si02 and 
Ta205 or (b) AI2O3 and Ta205. Shown are loss factors for frequencies up to 100 kHz. 

a value for the specific heat capacity of Ta205 (bulk) from No values for thermal 

conductivity or Poisson's ratio were available for Ta205; absent better information we take 
these to be closer to those of sapphire than silica, and assign them the same values as 
used for AI2O3. Using Eq. IHlthe thermoelastic losses from coatings of either Si02/Ta205 
or Al203/Ta205 applied to silica and sapphire substrates were calculated for frequencies up 
to 100 kHz, a typical range of interest for measurements of mechanical loss. The estimated 
loss factors are plotted in figures [T] and El 

It can be seen from these figures that in general, the calculated magnitude of the ther- 
moelastic losses from these mirror and coating combinations can be of the order of a few 
10^^ to approximately 10~^, comparable to the levels of coating loss factors predicted to be 
significant in estimations of the thermal noise level in advanced gravitational wave detec- 
tors J 3. 

In previous work measurements have been made of the mechanical loss in the fre- 
quency range from 2.8kHz to 73kHz, for dielectric coatings of Si02/Ta205 applied to fused 
silica substrates. The measured coating loss factors were found to be of the order of 2.8 x 
10~^. It can be seen from figure ^ curve (a) that the estimated thermoelastic losses for this 
particular coating/substrate combination are much smaller than the measured losses. This 
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2 4 6 8 10 

frequency (Hz) 

FIG. 2: Frequency dependence of the thermoelastic losses of sapphire substrates with dielectric 
multilayer coatings formed from thirty alternating quarter-lambda layers of (a) Si02 and Ta205 
or (b) AI2O3 and Ta205. Shown are loss factors for frequencies up to 100 kHz. 

suggests that the measured losses are not predominantly thermoelastic in origin, and are 
associated with some other form of dissipation. 

However, the thermoelastic losses for other combinations of mirror and coating mate- 
rials are estimated to be considerably larger than is the case for Si02/Ta205 coatings on 
silica substrates, see for example figure ^ curve (b) for an Al203/Ta205 coating on a silica 
substrate. Thus this form of dissipation should be considered in the interpretation of mea- 
surements of coating loss factors. Since this form of dissipation is frequency dependent, it 
is clearly important to estimate the magnitude of the dissipation in the frequency range of 
interest for gravitational wave detectors. Figure El plots the same loss factors shown in the 
figures above, focussing on frequencies up to approximately IkHz. 

From figureiniit can be seen that the thermoelastic loss in the gravitational- wave detection 
band is lower than at the higher frequencies sampled by mechanical loss measurements, 
however it can still be of the order of 10~^. It should be noted that in the absence of 
dissipation from other sources thermal noise arising from coating-related thermoelastic losses 
will form a limit to the thermal-noise performance of interferometric detectors in a manner 
similar to the thermoelastic noise from the substrates themselves, y Section ITVDl thus 
addresses the derivation of an expression for the thermal noise from coatings arising from 
thermoelastic dissipation, numerical results from which are presented in section IIII Bl 
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FIG. 3: Calculated frequency dependence of the thermoelastic losses of (a) an Si02/Ta205 coating 
applied to a silica substrate (b) an Al203/Ta205 applied to a silica substrate (c) an Si02/Ta205 
coating applied to a sapphire substrate and (d) an Al203/Ta2 05 coating applied to a silica sub- 
strate. Shown are the loss factors for frequencies up to IkHz. 

B. Thermal noise 



Using Eqs. ^1 and [721 with the parameters for coating and substrate properties given ear- 



lier, the thermal displacement noise resulting from thermoelastic dissipation, ySx{f), asso- 
ciated with silica and sapphire mirrors with coatings of either Si02/Ta205 or Al203/Ta205 
can be estimated. Here, multilayer coatings of a thickness equivalent to 10 ppm transmission 
were modelled, since this represents a typical specification for the transmission of a mirror 
coating used in the Fabry- Perot arm cavities of a gravitational wave detector. Figure E] 
shows the calculated noise for each case, for frequencies up to 1 kHz. A beam radius, w, of 
5.5 cm was assumed. 

For comparison, the target level for total displacement noise per test mass in the Advanced 
LIGO gravitational wave interferometer design is approximately 6 x 10~^^m/-\/Hz at 100 Hz 
if sapphire mirrors are used and approximately 8x10 ^^^m/ v^Hz for silica mirrors. Q Fig- 
ure m shows that for Si02/Ta205 coatings on sihca or sapphire substrates and Al203/Ta205 
coatings on sapphire substrates, the expected coating-related thermoelastic displacement 
noise is below the required specification at 100 Hz. However for an Al203/Ta205 coating 
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frequency [Hz] 



FIG. 4: Calculated thermoelastic thermal noise, \/Sx{f) for coatings of (a) an Si02/Ta205 coating 
applied to a silica substrate (b) an Al203/Ta205 applied to a silica substrate (c) an Si02/Ta205 
coating applied to a sapphire substrate and (d) an Al203/Ta2 05 coating applied to a silica sub- 
strate. Shown in each case is the noise obtained using the explicit low frequency limit for the noise 
(solid line), and the noise including the full frequency dependence (dashed line), as given in Eg. 1101 

on a fused silica substrate the noise from this dissipation mechanism alone is above the 
specification at 100 Hz. It is also clear that the same coating will result in a different level 
of noise if applied to different substrates. 

There are several other points illustrated by figure E] worth consideration. Firstly, for 
some coating/substrate combinations the thermoelastic noise starts to deviate significantly 
from the explicit low frequency limit in the frequency range of interest for detector operation. 
Thus when estimating the expected level of this noise it is important to use the full frequency- 
dependent expression. In addition, comparing figures El and El it can be seen that whilst the 
thermoelastic loss for strain fields associated with typical loss measurements, (pi^ \\ is higher for 
a Si02/Ta205 coating on a sapphire substrate than for that coating on a silica substrate, the 
opposite trend holds for the thermoelastic displacement noise sensed in an interferometer. 
This seeming contradiction can be understood by comparing Eqs. EZl for 0;^ y and ^2 for 
0z,_L, the loss for a specified surface-normal stress. We see that the dependence of these 
two loss coefficients on the material properties is quite different, in particular containing a 
ratio of Young's moduli in the former but not the latter. In fact, for the material properties 



14 



characteristic of these coatings, unhke follows the same trend as Sx, consistent 

with the observation that the axial stress is large in the region of high optical intensity for 
the fields of Eq. 1^21 used in the noise calculation. 

It is also important to note that in each case the exact level of thermoelastic noise 
is a strong function of certain of the material parameters of the coatings, in particular 
the coefficient of thermal expansion, and given the lack of information available on the 
thermoelastic properties of ion-beam-sputtered coatings, our calculations here were carried 
out using plausible rather than definitive values for relevant coating material parameters. 
Thus these figures should be taken as estimates of the expected thermoelastic noise due to 
the coatings, rather than reliable results that can be used in design calculations. 



IV. DETAILED CALCULATION 



A. Thermal fields 



We start by calculating the thermal field 9{z,t) generated by the applied zeroth-order 
elastic fields from Appendix ^ We can cast the one- dimensional heat equation in the 
form 

'dt " ""'^ ~ "(l-2z/,)C,-9t^^°'"'^' ^ ^ 

where 9j is the time- varying temperature, = kj/Cj is the thermal diffusivity, Ej is the 
Young's modulus, Vj is the Poisson ratio, T is the background temperature, Cj is the heat 
capacity per unit volume, eo,M,j is the zeroth-order z-polarized compressional strain, and 
j = f,s indicates quantities evaluated in the film and the substrate, respectively. Taking 
sinusoidally time varying quantities according to a{z, t) = Re[a{z) exp{iut)], Eq.lT^ becomes 

iu9j{z) - Kj = -iujpj (14) 

where the source term is proportional to 

The quantity representing the sum of the strains according to 
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= ^ ^o,Mj ) (16) 



i=l 
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is proportional to the zeroth-order elastic field's amplitude with a combination of elastic 
constants that depends on the specific case, and is evaluated in Appendix |X1 a/ is possibly 
z-dependent to allow for spatially varying thermal expansion coefficient within the film. 
The boundary conditions are zero heat flux aX z = 0, continuity of heat flux at z = I, and 
vanishing heat flux for 2; — oo, or 



dz 



2=0 



dz 



z=l 



dz 



de^ 

dz 



0, 



(17) 



respectively. 

The total solution can be constructed as the sum of homogeneous and particular solutions. 
Homogeneous solutions meeting the boundary conditions at 2; = and z ^ 00 are 



9hj{z) = 9if cosh(7/z) and 9h,s{z) = due 



(18) 



where the complex propagation constants of the damped thermal waves in the film and 
substrate are 

' (19) 



Ij = (1 + 0\/^/2/ti 



and 6if and 6is are constants determined by the boundary condition at z = I and the 
particular solution. 

The particular solutions will be evaluated for two specific cases of practical interest in 
Appendix ini For the time being take them to be 9pj{z), so that the total solutions are 



^/(^) = ^P,/(^) + ^1/ cosh(7/ z) 



(20) 



Note that both 0pj{z) and 61 j will be proportional to the amplitude of the zeroth-order 
elastic fields. 



B. Elastic fields and energy density 

The rate at which work is done per unit volume by internal stresses on a deformed body 



IS 



P = a.,^ . (21) 



This expression is correct independent of whether the body responds elastically or anelasti- 
cally to the stresses. The cycle average of the delivered power density (or, equivalently, 
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the average dissipated power density) is then, for fields sinusoidal in time of the form 

aiij{z,t) = Re{[ao,ii,j + cri^ujiz)] exp{iujt)} and similar for eiij{z,t), 

u ^ 

Vdiss{z) = -lTl]lm[(T*;.£ife] 
^ i=l 
^3 

^ i=l 

^ r 

~ -irlllm [(T* iiei,,i + ai*iieo J . (22) 
^ i=i 

The last form of this equation is justified by the following observations: For our problem, only 
the longitudinal strains are significant and the zeroth-order elastic fields are real (Appendix 
1X1 and Eqs. while the first-order elastic fields (those that depend on the thermal 

field. Appendix [nj are complex and so contribute to the dissipation. We assume that the 
dissipation is small, so that a calculation to lowest order in the thermal field will be adequate, 
and second order terms can be dropped. 

To evaluate Eq. |221 we need the zeroth- and first-order elastic fields in the film and the 
substrate, as derived in Appendices El and O It is convenient to write the zeroth-order 
fields in terms of Uo and two combinations of the fields that are invariant through the 
region of interest in the body. These are defined by (Tq = o"o,22, the compressional stress 
normal to the surface of the object, and = {.^o,xx + £^o,j/s/)/2, the in-plane dilation. For 
convenience, we can take Eq^xx = ^o.yy, though only their sum matters for the thermoelastic 
calculation. The antisymmetric combination of in-plane strains is a pure shear, does not 
interact thermoelastically, and so be can neglected in this analysis (as is discussed in section 
live Note that in cylindrical coordinates Eq = (e^rr + £^o,6»e)/2. The zeroth-order fields 
can then be summarized as 

where j = /, s, and the (^o,M,i and -Bo,ii,j) and the {ao^aj and bo n j) are combinations of 
elastic constants given in Eqs. I A 71 and IA141 respectively. Similarly, it is convenient to write 
the first-order elastic fields as proportional to the local temperature and thermal expansion 
coefficient, since we assume that the frequencies of interest are low enough that the elastic 
response can be treated quasistatically: 

Ei^iij{z) = Ai^iijajejiz), ai^iij{z) = B^^.i^^ajOjiz) (24) 
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The (real) coefficients Ai^uj and are given in Eqs. IC7I 

With Eqs. 1^ and 1^ in Eq. |221 the dissipated power density can be written as 




(25) 



where 



3 



= H [Bo,ii,jA 




i=l 



3 




(26) 



i=l 



The dissipated power per unit area is then obtained by integrating the power density over 
the thickness of the body, 



where we allow for the possibility of spatial variation in the thermal expansion coefficient in 
the film but assume that it is uniform in the substrate. 

We can use this expression both to evaluate the dissipation that would be measured in 
a typical measurement of the Q of a coated test mass and to calculate the coating-related 
thermoelastic contribution to the displacement noise imposed on an optical field incident on 
a test mass. The most convenient form of the analysis is somewhat different in these two 
contexts. We begin with the case of a Q measurement. 

C. Effective thermoelastic losses in measurements of elastic Q 

While the total dissipated power given in Eq. 1221 is the physical quantity of importance to 
measurements of Q, and is nonlocal in nature, occurring in both the film and the substrate, 
it generally occurs in a region thin compared to the dimensions of the test mass, so for 
comparison with experimental results it is convenient to describe the loss in terms of an 
effective associated with the coating. To define such an effective 0, we must compare 
the dissipated power to some stored energy. A reasonable choice of stored energy for the 




(27) 
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definition of d) is that in the film, i.e. 



I 3 

[/stor/area = ttXI^^ cro,n,/^o,ii,j 



2 . 

1=1 



V-fUf, (28) 



where 

3 



f^/ = E^o..,Mo..,/. (29) 

i=l 

In writing Eqs. |2H1 and EHl we assumed that the film is on a stress-free surface, so that oq^^z 
and hence ctq vanish. Uf is calculated in Appendix [XJ Eq. IA9I We then have for 0, with 
Eq. 1211 

-Pdiss^ 

2vrf/stor 

= 0/ + 0« (30) 



where 



and 



^ faf{z)lm[-9f{z)/eo]dz (31) 

U fL Jo 



^^-^^^as lm[-es{z)/eo]dz . (32) 
U fl Jl 

Note that since the thermal fields are proportional to Eq, the quantity in square brackets in 
Eqs. EH and 1221 is independent of Eq, as are all the other factors in these equations. 

To make further progress, we must find the particular solution and the coefficients for 
the homogeneous solutions for the thermal field in the specific cases of interest. We consider 
the specific cases of a uniform film on a uniform substrate, and a periodic film on a uniform 
substrate. In Appendix El we obtain the thermal fields for these two cases. 



1. Uniform film and substrate 

Consider first the simple model of a uniform film on a substrate, with possibly different 
thermophysical properties in film and substrate. For this case, we can take the thermal 
expansion coefficients in the film and substrate to he af{z) = af and a^, respectively, and 
the particular solutions to the heat equation, given in Appendix |Bl as Eq. IB5t are 

e,,{z) = -(3,. (33) 
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Since the particular solutions are real, the only contribution to the imaginary part of the 
integrals in Eqs. |^ and 1221 come from the homogeneous solutions, so with Eqs. ^Iwe have 
for the film, 



D 



f 



Ufl Jo 
Dfaj 



af{z) lm[—9f{z)/eo] dz 



Ufl 
Df(^f 



Im 
Im 



-(^ij/^o) / cosh(7/z)c?z 



-(^i,//£o)7/^sinh(7j/) 



and similarly for the substrate 



Im 



(34) 



(35) 



where the Dj are combinations of elastic constants defined in Eq. |26l and calculated in Eq. 
IC8I Summing the contributions to the dissipation from the film and the substrate, Eqs 
and we can express the total loss as 



Ufl 



Im 



-(^i,//£o)7/^sinh(7j/) 



H Im 

Ufl 



(36) 



where the subscript / is used to indicate a quantity resulting from thermoelastic behavior 
over the entire thickness of the film (contrasted to multilayer case in following section) and 
II indicates the case of specified in-plane strain (in contrast to _L for specified surface- normal 
stress). With the coefficients 9ij from Appendix FBI Eqs. IB6I and IE7I inserted into Eq. inSl 
we can express the total loss, after some algebra, as 



eo Ef 
2Cf{l - Uf)T 

2Efa}T 



1 - Uf 

afEf 



Vs Is 

a^E, 



1 2 



1 2 



a,Es{l-Uf)Cf 



Cf{l-Uf) [ afEf{l-u,)C,\ 
where the frequency dependence is contained in the function g{uj) defined by 



(37) 



Im 



smhL iujTf) 



(38) 



In deriving this result, we made use of Dj from Eq. ICSl Uf from Eq. IA91 R from Eq. IB41 
A/5 = Pf — Ps from Eq. [T31 S from Eqs. lAll and lA81 and defined r/ = I'^/nf so that with 
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Eq. ^Jwe have 7// = ^iuTj. Note that the combination of material properties in square 
brackets in Eq. EZlis positive definite and vanishes if the film and substrate properties are 
identical. The quadratic dependence on the difference between substrate and film properties 
can lead to dissipation that is sensitive to small changes in the film properties. 

A similar analysis can be carried out for an antisymmetric in-plane strain, So^xx = —£o,yy 
We find that = 0, so that no thermal wave is generated (consistent with the observation 
that this antisymmetric strain is a pure shear along axes rotated 7r/4 with respect to x 
and y, causes no volume change, and hence does not contribute to the source term for the 
thermal wave). We also find that Dj = 0, indicating that there will be no power dissipated 
by the interaction of the zeroth-order antisymmetric strain with the first-order strain fields 
generated by the thermal wave (driven by other zeroth-order strains possibly present). This 
latter observation can be explained by noting that the thermal wave generates no first-order 
shear strains (for the geometry considered here), and that there is no energy term associated 
with the product of shear and compressional strains in isotropic media. 

While not encountered in the measurement of elastic loss in coated masses, an expression 
for the dissipation for a specified surface-normal stress is useful for developing an understand- 
ing of the thermal noise results in comparison to results for loss measurements. Following 
the same analysis as was used to find (pi 11, but replacing the stored energy in Eq. OHl with 



t/stor/area = l—uj 



(39) 



where 

we find that replacing Uf 

for 0; I , 



(40) 



1=1 



Uf, Dj dj, and Eq in Eq. 123 yields the correct result 



1 + 



as (l-z//)(l + i/,)C/ 



1 2 



(41) 



Cf (1 - - 2uf) [ af{l- z/,)(l + Uf) Cs_ 

where we made use of dj from Eq. I(^9[ Uf from Eq. IA16[ E from Eqs. lAll and IA15| other 
quantities as after Eq. EHl and g{uj) is the same frequency dependence given in Eq. 
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2. Modulated film and uniform substrate 



Optical coatings of interest for use on test masses for gravitational wave detectors are 
invariably multilayered, so the analysis in the previous section of a uniform film cannot 
be correct in detail. Because the thermal diffusion length for frequencies of interest are in 
general long compared to the thickness of individual layers in the film, an accurate approx- 
imation for a multilayer coating can be obtained by using a suitable averaging process to 
model it as a uniform layer (except at very high frequencies), as discussed in Appendix iDl 
Following the procedure described there, we find that the result in Eq. ^7\is replaced by 

2 

g{uj) (42) 



2CfT 



1 / Ea \ 1 E^a^ 



Ci? \1 — Z//avg C*s 1 



and the frequency dependence g^uj) is unchanged except for replacing the time constant r/ 
by an appropriately averaged one tf- The volume- weighted average indicated by (Ar)avg is 
defined in Eq. ID 11 the averaged heat capacity Cp in Eq. ID41 and Tp in Eq. ID61 

While the results based on this averaging process appear reasonable, it is useful to explore 
for unexpected effects associated with the spatial variation of thermoelastic properties within 
the multilayer film. As a simple model of such a case, we consider a film whose thermal 
expansion coefficient, but no other property, varies periodically, on a uniform substrate. The 
calculation is similar in principal, but somewhat more complicated than for the case of a 
uniform film on a uniform substrate. For this case, we take a thermal expansion coefficient 
of the form 

af{z) = af + amCos{Kmz) and as{z) = as (43) 

and the particular solution to the heat equation in the film, given in Appendix ^ as Eq. 
iBTal is 

^P,/(^) = -Pf - Pm 2 1V2 C0S(JS:„Z) (44) 
7/ + 



while the particular solution in the substrate, Eq. IB 141 remains the same as for the uniform 
film case 

Op,s = -Ps . (45) 

For simplicity, we consider here only the case of specified in-plane strain and vanishing 
surface- normal stress ctq. Since the particular solution in the film has a complex part, the 
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imaginary part of the integrals in Eqs. |M] and |221 will contain contributions from both the 
particular and the homogeneous solutions. With Eqs. ^1 021 and El we have in the film 



where 



Ufl 



fl JO 



Oif{z) lm[—6f{z)/eQ\ dz 



f,h) , 



- / [af + amCos{Kmz)] [9pj{z)/eo] dz 



m I 



Km ^0 



Kml ^ sm{2Kml) 



Im 



7| 



(46) 



(47) 



and 



^f,h = Im 



— Im 



+ a, 



- / [af + aniCos{Kmz)] [6h,f{z)/eo] dz 

•J 

af — -'^f smh7// 
6i f Km cosh(7j/) sin(_ft'm/) + 7j sinh(7j/) cos(_K'm/) 



(48) 



£0 + 7/ 

Since the thermal expansion coefficient in the substrate is assumed uniform, the result 
for (ps is similar to that of Eq. jSHl (though of course ^i^s will be different in the two cases). 
We use Eqs. ITHl EHl and E31 to obtain for the substrate 



Ufl 



Im 



-(^i,s/£o)7s e 



(49) 



Considerably more effort is required to convert Eqs. 03 EHl and EHl into a simple form like 
Eq. EZl after inserting 9ij and from Eqs. IBlSl and IB16I The following two terms emerge 
as dominant in Eqs. 1171 and ilHl for cases where \Km\ » |7/|, 



so that with Eq. 



we have 

Df J ami f3„ 



O^ml f-^m, J 

af—g{uj) , 

^0 



Im 



7| 



Ufl [ 2 so [jj + Kl 
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+ «/ — gM 

^0 



(50) 



(51) 



where g{uj) is defined in Eq. EH Similarly, keeping the leading term in Eq. 0^ after inserting 
Eq. lBT6l leads to 

-Rg{uj) (52) 



£q Ufl 

where all neglected terms are smaller by at least one factor of '-/f/Km than those retained. 
The first term in Eq. |^ is unique to a modulated film. The second term in Eq. |^ and 
Eq. 1221 are just the same as those that appeared in the expression for a homogeneous film 
and substrate, Eq. EHl with «/ — a/, so that portion of the solution can be used here 
immediately. Rewriting Eqs. |^ and with Dj from Eq. IC81 Uf from Eq. IA91 R from Eq. 
IM I3m from Eq. lETOl I3j from Eq. US 7j from Eq. d and from Eq. HUwith Eq. \M 
we obtain 



7/ 



+ 0Mh (53) 



where 0/, \\ is the loss for a uniform film on a uniform substrate given in Eq. EZl with 
af Of. Defining a characteristic time for diffusion in the modulated structure, 

"""^ " ^ ^ ^ ^^^^ 
where the second form follows from the definition of r/ following Eq. EH 4'tot, \\ can be written 
in the form 



where the frequency dependence gm{'^) is 



g^{uj) + (pi^W , (55) 



3m{^) = 1 ^'^"l 2 ■ (56) 
1 + uj^ri 

In section llV C 31 we show that for the frequencies of interest, the frequency dependence 
(7m (u;) from Eq. IKHl representing heat flow between the multilayers is, as expected, small 
compared to g{uj) from Eq. IHHl representing heat flow between film and substrate, so that 
considering only the contribution of the averaged form 0;^ y from Eq. 021 is a good approxi- 
mation. 

3. Frequency dependence 

The functions g{uj) and gm{^), defined in Eqs. EHland|Sni respectively, contain all the 
frequency dependence of the dissipation, and will reappear in our discussion of displacement 
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noise in section ITVDl It is therefore worth investigating their general features in some detaiL 
Consider first g{u!)] it is convenient to define the real variable ^ according to 



7// = y^=(l+0e/2 



(57) 



so that 



2uJTf 



(58) 



where from Eq. |2H1t/ = ^"^Cf/kf . The frequency dependence of the loss due to a uniform 
film can then be written as 



9(0 = 



sinh^ — sin^ + _R(cosh^ — cos^) 
cosh^ + cos^ + 2i?sinh^ + i?2(cosh^ — cos^) 



(59) 



R = 
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- - R = 
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10 




FIG. 5: Frequency dependent part of thermoelastic loss function, g{^) as a function of frequency, 
where ^ = y^TurrJ and rj = PCf/kj. Curves are shown for three values of the parameter R = 



kfjf/kg'js = y kfCf /kgCs- Note that peak height is only a weak function of R. 

In terms of the normalized frequency ^, the normalized dissipation depends only on one 



parameter, R, defined after Eq. IB4I as i? = kf'jf/ks'js = ykfCj/ksCs- Figure El shows the 
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0.45 



0.35 




0.25 - 



FIG. 6: Peak value of the normalized dissipation, gmax = ^(Cmax) as a function of the parameter 
R. The peak value of the normalized dissipation is seen to be only a weak function of R, the only 
material parameter on which it depends. 



frequency dependence of the dissipation for R = (0.1, 1, 10). Figure El shows the dependence 
on R of gmax, the peak value of the normalized dissipation, and figured the dependence on 
R of ^max5 the normalized frequency at which this peak occurs. We see that g^^.^ depends 
only weakly on R, ranging from 0.41 for i? <^ 1 to 0.21 for i? ^ 1. ^ma^ is close to the 
thermal diffusion time across the film = 2) for R < 1, and decreases as 1/i? with R for 
R>1. 

Useful forms for low- and high-frequency limits of the dissipation can be obtained from 
Eq. EHl Expanding for ,^ <^ 1 we find 



9(0 



while for ^ ^ 1 



9(0 



(60) 



(61) 



Since u oc the leading behavior for low frequencies goes as the i/tJ, and, surprisingly, the 
sign of the term linear in frequency depends on the value of R, crossing zero for R = l/VS. 
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FIG. 7: Normalized frequency ^max at which the normaHzed dissipation takes its maximum value, 
as a function of the parameter R. ^max corresponds to the thermal diffusion time across the film 
for < 1, and decreases as 1/R (corresponding to a;inax decreasing as l/R"^ for R> 1). 

At high frequencies, the dissipation falls off as 1/ ^/uJ. 

The frequency dependence described by gm{'^) for the contribution from the multilayer 
coating is simple, and essentially the same as that for conventional thermoelastic damping, so 
it needs little further discussion. It is important to note that for typical multilayer coatings, 
the characteristic time for the multilayer effects is much shorter than t/ for the effects of 
the averaged uniform layer; for a coating with 2N layers, 

Tj/Tra = K^f = Att'N' . (62) 

Since a typical high reflector that might be used in a LIGO interferometer has 40 layers, t/ ~ 
IGOOOTm, so that the peak frequency for the contribution to the thermoelastic dissipation 
from thermal diffusion between the layers will be at a frequency ~ 16000 times higher, 
generally pushing the peak well above typical measurement ranges. The thermoelastic effects 
at frequencies of interest either for elastic Q measurements or thermal noise are thus generally 
dominated by the contributions of the averaged film. 
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D. Thermal noise 



The results of section UlIAI indicate that the power dissipated by the thermoelastic effects 
can be comparable or even exceed that dissipated by the elastic loss in typical multilayer 
coatings. It is then reasonable to assume that the magnitude of the noise induced by the 
thermoelastic mechanism could be comparable to that from the elastic loss, and therefore 
must be calculated as part of the total noise budget for the coated mass. 

Following the approach of Q], the displacement noise imposed on a Gaussian beam of 
normalized intensity distribution /(r) 

/(r) = ^expf^) (63) 



SAf)df= J,, S' df (64) 



is given by 

2kBT Wdis s 

where Wdiss is the cycle-averaged power dissipated by a pressure field p(r) oscillating at a 
frequency uj = 27r/, of the same radial distribution as the intensity and with a resultant 
force Fq, i.e. 

p(r) = FqHt) cos(cut) . (65) 

Since the radius of the Gaussian beam is much larger than the thickness of the film or the 
thermal-wave decay length, we can use the one- dimensional theory developed in previous 
sections of this paper to evaluate PVdiss- The zeroth-order elastic fields required for this 
calculation are available in j^. We have from their Eqs. (AlO) 

= ~pir) 

5o(r) = [err{r) + e60{r)]/2 

= -p(r)/4(A + /i) 
(l + i.,)(l-2z/,) 



p(r) (66) 



2Es 

where the third form of follows from the definition of the Lame constants in terms of the 
Young's modulus and Poisson ratio. 

Here we analyze the noise due to an appropriately averaged uniform film, since the 
analysis of section ITV G 31 showed that the contribution to the dissipation associated with the 
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thermal diffusion between the layers within the film are significant only at frequencies well 
above the LIGO detection band. We first consider a uniform film, and then the modifications 
necessary to describe an appropriately averaged multilayer. 

Starting with Eq. [27|for the dissipated power per unit area, and Eqs. 1201 foi' the thermal 
fields we have 



diss 



area 



= - {afiDfEo + dfao)lm[-9ifjf ^ sinh(7//)] 

(-7/0]} 



ujA(3l 



Cf 



(67) 



where as in Eqs. |^ and 1221 we simplified the result by noting that only the homogeneous 
part of the thermal solutions contains an imaginary part, and in the second form used Eqs. 
IB6I and IB 71 for 6ij and ^i .,. The frequency dependence is contained in the same function 
g{uj) defined in Eq. EHl 

For this case, where there are both axial stresses and in-plane strains, we can calculate 
Ap = pf- (3s from Eqs. US CHI HB |Ml andETB We find that 



C, 1- 
Ci 1 - 



1 + z/,- 
2^0 + —ET^^^o 
Ej 

■(l + z/,)(l-2z/,) , l + Vj 



+ 



Eg Ej 



(68) 



where the second follows for the specific form of the elastic fields given in Eqs. EHl With 
Eqs. IC8l and [C9l for the combinations of elastic constants represented by Dj and dj, and A(3 
calculated from Eq. IHHl the result in Eq. IHZIfor the dissipated power per unit area becomes 
Pdiss('^) _ ujTlCf 



area 



X 



af 



-p(r) g{uj) 



- ^2(1 + Us] 



i-.r 1-., Es 

where we replaced cxg with p(r)^ according to Eq. ESI Integrating over the infinite cross- 
section to obtain the total dissipated power VFdiss, and inserting that result into Eq. EH for 
Sx{f), we finally obtain 



X 



aj 1 
.^2 



Es 



[l + l^s 



(70) 
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where is a dimensionless positive-definite combination of material constants tliat vanisli 
wfien tlie film and substrate are identical, 

^ 2 



2asCf{l-Vf) 



l + Vf 



1 



(71) 



Eq. [701 is the final result for the thermoelastic displacement noise associated with a uniform 
coating. The frequency dependence represented by g{uj) is discussed at length in section 

An accurate approximation for a multilayer coating can be obtained by using a suitable 
averaging process to model it as a uniform layer, as discussed in Appendix |D| Following the 
procedure described there, the result in Eq. [TUlis replaced by 



where 



A^ 



2a,Cf 



7r2/ ^ C2 ' 
l + v 



1 + Vsf~^^9{uj) 



a 



+ 1 - 



E 



1 



(72) 



(73) 



avg 



and the frequency dependence giu) is unchanged except for replacing the time constant 
by an appropriately averaged one Tp. The volume- weighted average indicated by (X)avg is 
defined in Eq. ID 11 the averaged heat capacity in Eq. ID41 and Tp in Eq. ID6I Since for 
room-temperature operation the thermoelastic noise is most important at frequencies falling 
in the low-frequency limit of g{uj), it is useful to insert into Eq. [721 the approximate result 
for g given in Eq. [Hllto obtain 

.,2 



n I J- -|- Z^s 



a: 



(74) 



Note that, as seen in figure HI this low frequency limit becomes inaccurate at the upper end 
of the gravitational-wave detection band. 



V. SUMMARY AND CONCLUSIONS 



In this paper we have derived expressions for the thermoelastic dissipation associated 
with a coating on a test mass. For strains of the type consistent with mechanical loss 
measurements, numerical evaluation of the thermoelastic loss factors for coating/test mass 
material combinations of the type being considered for use in future gravitational wave in- 
terferometers shows that thermoelastic dissipation is of a level comparable to that predicted 
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to affect the sensitivity of advanced interferometers. Also derived is an expression for the 
expected power spectral density of thermoelastic noise from the coating of a mirror inter- 
rogated with a Gaussian beam. Evaluating this expression across the gravitational wave 
detector band using plausible values for the material parameters of coatings and substrates 
results in displacement noise that in some cases exceeds typical design sensitivities. 

It should be noted that the expected thermoelastic noise is a strong function of the 
difference of the material parameters in the substrate and coating, so that the same coating 
will have different thermoelastic losses on different substrates. As many of the necessary 
material parameters are not well characterized, the noise levels calculated here should be 
considered as estimates only. Further experimental measurements of coating dissipation 
for likely choices for coating and substrate materials, and better characterization of the 
intrinsic coating thermophysical properties, should allow more accurate determination of 
the magnitude of the thermoelastic effects. 
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APPENDIX A: ZEROTH-ORDER DRIVING FIELDS 

We need a set of zeroth order driving fields that are consistent with the elastic boundary 
conditions. In all cases we must have continuity of the in-plane strains and the normal 
stress at the film-substrate interface, z = I: eo,xx,s = £o,xxj, £o,yy,s = ^o,yy,f and cro,zzj = 
o'o,zz,s- As discussed in section HI the pertinent elastic fields can be specified in terms of 
two components, the in-plane dilation Eq = {eq^xx + £o,yy)/'^ and the axial stress ctq = ctq^zz, 
which are independent of z under the assumptions set up in section |lj We can neglect the 
anti-symmetric in-plane strain {eq^xx — ^o,yy) which does not interact thermoelastically (as is 
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shown in section ITVC and, for convenience, can take Eq^xx = £o,yy = £o- 

We consider two cases, a stress-free surface with a specified in-plane strain (ctq = and 
£0 specified), and a specified surface- normal stress with a vanishing in-plane strain {eq = 
and (To specified). Any elastic state pertinent to the thermoelastic problem can be obtained 
as an appropriately weighted sum of these two solutions. For the general case, where both 
(To and Eo are nonzero, one of the important results of this appendix, Ej defined in Eq. 
can be written 

3 

— ^ ] ^0,ii,j 
1=1 

= SjEo + sjao . (Al) 

The combination of elastic constants Sj and Sj are obtained in this Appendix, Eqs. IA8I and 
IA15| respectively. 



1. Specified in-plane strain, stress- free surface 

For a stress-free surface of the mass, as would be the case for a Q-measurement, we have 
(To = and Eq specified. Noting that the under these assumptions the continuity condition on 
the normal stress implies that crQ^zz,s = <^o,zzj = ctq = 0, the only unknown field components 
are ^o.zz,/; ^o,zz,s; (^o,xxj = '^o,yyj = '^0,11,/) and <Jo^\\^s- The symmetry of the problem allowed 
us to take (To,xx = cro,3/y = CTo.ii- 

Begin with Hooke's law, Eq. 5.14 of |l8j : 

£0,zz,j = ■^ Wo,zz,j - ^j{(^0,xx,j + Cro,yy,j)] 

= -^^0,11. • (A2) 

where we recall the notation that a subscript j = f,s stands for a quantity evaluated in the 
film or substrate, respectively. Summing the expressions for the in-plane strains in it 
follows that 

£0,xx,j + £0,yy,j = -TrWo,xx,j + (^0,yy,j " ^j{(^0,xx,j + (^0,yy,j)] (A3) 
Ej 



or equivalently 



^0 = —Er^(^o,\\,j ■ (A4) 
Ej 
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Finally, going back to Eq. IA2I with Eq. IA41 we find 



^0,zz,j — —^0 ■ 

It is convenient to summarize these results for the zeroth-order elastic fields in the form 
given in Eq. 1221 

^0,Mj = ^0,Mj^O; '^O.iij = Bq^hjEo, (A6) 

where 

A -A - 1 4 - "^'^J 

■^(},xx,j 0,yy,j ^1 ^0,zz,i -, 

i — Vj 

E ■ 

i — Uj 

A result, used in Eq. UHl is the evaluation of a sum over strains introduced in Eq. ITHj 
which, with Eqs. lAll and I A 71 becomes 

3 3 2(1 — 2z/) 

i=l i=l ^ " ^3 

Another result, used in Eq. 1211 to evaluate the energy stored in the film is 

3 2Ef 
Uf^j:Bo,ujAo,uj = -^ . (A9) 
i=i ^-^f 

2. Specified surface-normal stress, vanishing in-plane strain 



For a specified surface-normal stress, as would be the case for calculating thermal noise, 
we have cro,zzj = o"o- To make this case complementary to that in section IXTl we assume 
vanishing in-plane strains, i.e. Eq^xx = £o,yy = eo = in both the film and the substrate. 
Noting that the under these assumptions the continuity condition on the normal stress 
imphes that o"o,^^,s = O"o,z2,/ = ctq, so the only unknown field components are So,zz,j and 
(Tq^xxj = <^o,yy,j = J = /) -S- The analysis is similar to that in section \KT\ Begin 

with Eq. 5.14 of [fj: 



^0,xx,j - ■Tr[(^0,xx,j - ^j{<^0,yy,j + (^0,zz,j)] (AlO) 

which can be solved with cro,xx,j = <^o,yy,j = (^o,\\,j to yield 

^o,||,i = T^fxo . (All) 

1 - Uj 
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With another of Eqs. 5.14 from 



£o,zz,j - -TrWo,zz,j - i^j{cro,xx,j + o"o,j/s/j)] (A12) 



and Eq. lAlll we obtain 

ao (l-2z/,)(l + t/,) 
Ej 1 - z/j 

We can again collect the results of this section in the form given in Eqs 



eo,zz, = . (A13) 



a-0,xx,j — O'O^yy,] — 0, O,0,zz,j i 771 

1 - Uj hj 

bo,xx,j = bo^yyj = - — - — , bo^zzj = 1 • (^14) 

i — Uj 

A result, used in Eq. UHl is the evaluation of a sum over strains introduced in Eq. [121 
which, with Eqs. lAll and IA14I can be written 

_A _A _ (1 - 2z/j)(l + z/j) 1 

i=l i=l 1 - ^1 

Another result, used in Eq. EUJto evaluate the energy stored in the film is 

x^, (1 - 2z/f)(l + z/f) , 

= E bo,^^Jao,^J = ^ Efil-Uf) ' ^^^^^ 

APPENDIX B: THE THERMAL FIELDS FOR TWO IMPORTANT CASES 

The unknown coefficients in the homogeneous parts of the thermal fields, Eqs. ^1 can 
be obtained from the particular solutions, Eqs. EOJ and the boundary conditions, Eqs. [T7| 
Continuity of the thermal field at 2; = / requires 

Opjil) + ^1/ cosh(7/ /) = ep,s + Oue-'^' (Bl) 

while continuity of the thermal fiux requires 



^;,/(0 + OiHf sinh(7y /) = -h eislsc-^'. (B2) 



Simultaneous solution of these equations yields 



^'^ cosh(7//) + i?sinh(7//) ^ ^ 



34 



and 

, ,j K,s - e,j{l)]RsmHjfl) + {Rhf) cosh(7//)g;,//) 

cosh(7y/) +i?sinh(7j/) ■ ^ ^ 

where R = kf'jf/ks'^s- To make further progress, it is necessary to find the particular 

solutions for specific cases. We consider here two cases of interest, a uniform film on a 

uniform substrate, and a periodic film on a uniform substrate. 

1. Uniform film on uniform substrate 

Consider first both the film and substrate to be uniform. By inspection of the thermal 
field equation, Eq. [T31 particular solutions for this case are constant and given by 

e,,{z) = -(3, . (B5) 

With Eq. IB 5 1 for the particular solutions, the coefficients in the homogeneous solutions in 
film and substrate from Eqs. IB3I and IB4I become 



''^ cosh(7y/) + i?sinh(7y/) ^ ^ 



and 



where A(3 = (3f - (3s 



= r^si Al3Rsmh{-ffl) 

cosh(7//) +i?sinh(7//) ■ ^ ' 



2. Modulated film on uniform substrate 

In the case of a nonuniform film, the expression for the particular solution is somewhat 
more complicated. Assume a film whose thermal conductivity takes the form 

af{z) = af + am cos{Kjnz) , (B8) 

in which case with Eq. takes the form 

(3fiz)=(3f + (3^cosiKmz) (B9) 

where 
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With Eq. IB91 the thermal field equation Eq. Uptakes the form 

= AfldpA^) +Pf + P^ cos{Kmz)] (Bll) 

where the definition 7^ = iuj/nf from Eq. ITTH was used to obtain the second form. The 
particular solution has two terms, a constant part similar to that in Eq. IBSI for the uniform 
case, and one that has a spatial variation that follows the thermal expansion coefficient. To 
obtain the periodic part, take an ansatz 9pj{z) = q cos{Kmz). With this ansatz in Eq. IBllI 
we obtain 

(7j + ir^)g = -7?/5n. (B12) 
Combining with the constant part we obtain the total particular solution in the film, 

9pj{z) = -(3f - ^m^P^ cos{Kmz). (B13) 

7/ + -f^-m 



For the assumed uniform substrate, the particular solution is like that in Eq. IB5t i.e. 

0,,s = -Ps . (B14) 
With Eqs. IB 131 and IB14I for the particular solutions, the coefficients in the homogeneous 



solutions in film and substrate from Eqs. IB3I and lB4l become 

n ^ ^ ,0 7/ cos(i^^/) - {RKmhf) smjKJ) 

''^ cosh{^ f I) + R sinhi^f I) '""^j + Kl cosh(7y/) +i?sinh(7jO ^ ^ 



and 

^1,, = - A(3e 



-yj i?sinh(7//) 

cosh(7/Z) + i?sinh(7/Z) 



7/ cosjK^l) sinh(7//) - {Km/if) sin(i^^/) cosh(7//) 



7/ + Kl^ cosh(7//) + R sinh(7//) 

where Ap = pf - p,. 



(B16) 



APPENDIX C: SOLVING FOR THE THERMOELASTICALLY GENERATED 
ELASTIC FIELDS 

Given the solution Eqs. EOlfor the oscillatory thermal field, we must solve for the ther- 
mally driven elastic fields, cri{z) and ei{z), whose imaginary parts lead to the dissipation 
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in which we are interested. The boundary conditions are ai^zz = at the stress-free surface 
2; = 0, si^xx = £i,yy = for z ^ 00, and continuity of the in-plane strains ei^xx and ei^yy, and 
the normal stress (Ji^zz, at the boundary between the film and substrate. 

The point of departure is the equation of elastic equilibrium, Eq. 7.8 of jlS^ . 

^ [Sxx + Syy + 2(1 - v)ezz - 2(1 + v)ae] = (CI) 

adapted here by dividing Landau's a by 3 to convert it from volumetric to linear expan- 
sion, replacing aW9 with V(a6') to accommodate a possible spatial variation in the thermal 
expansion coefficient, and specializing to stress and strain fields that depend only on z. 
Hooke's law in the presence of a nonuniform temperature field 9{z)^ Eq. 6.2 of 18j is: 

(C2) 



E 



V 

1 - 2z/ 



1 + z/ 

With the boundary condition a\^zz = 0, Eq. 1(^21 results in 



2u 



£l,zz,j - --, _\ i^hxxj + £l,yy,j) + i ,/ ^^j^j ■ (^3) 

Inserting Eq. K^SI into Eq. K^H we find 

Noting the continuity of the in-plane strains, and their vanishing at infinity, we conclude 
{.^i,xx,j + ^i,yy,j) = 0. With this result, Eq. IC3I becomes 

ei,zz,j = y^c^jOj . (C4) 

With the Hooke's law expression for axx analogous to Eq. IC2I for azz, and inserting Eq. IC4|, 
we obtain 

<7i,xx,j = — TZTT" (C5) 

and by symmetry ai^yyj = ai^xxj- 

These results constitute a consistent set of first-order elastic fields. It is convenient to 
summarize them in the form: 

= Ai^iijajOj, ai^iij{z) = Bi^iijajOj, (C6) 

where 



-'^i,xx,j ^,yyj ' ^ 



^l,xx,j — ^l,yy,j — 1 _ ,. 
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, B,^zz,j = . (C7) 



Combinations of these parameters used in calculating the dissipated power, Dj and dj in 
Eqs. |2ni can be evaluated with Eqs. I A 71 and IC7I as 

3 2^ . 

= [Bo,ii,jAi^iij - Bi^iijAo^iij] = ——^ , (C8) 

and with Eqs. IA14I and lCTI as 

3 1 + z/- 

dj = [bo,ii,jAi^iij - -Bi,njao,ijj] = _ ^ ■ (C9) 
i=i '^j 

APPENDIX D: AVERAGING MATERIAL PROPERTIES IN A PERIODIC 
MULTILAYER 

In cases of practical interest, the optical coating is a multilayer rather than a homogeneous 
film. The analysis in sections llV C 21 and II V C3l indicated that for realistic cases either for 
Q measurements or for thermoelastic noise, the thermal diffusion length is large compared 
to the period of the multilayer. Hence, an analysis that treats the film as an effective 
homogeneous medium with suitably averaged properties should yield a result of adequate 
accuracy. It is then necessary to form the appropriate average of the various material 
properties involved. For simplicity, we take the coating to consist of alternating layers of 
two types of material, labelled a and 6, of thicknesses da and db-, respectively. We define the 
volume averaging operator by 

(^X)^^ = Xg + Xb (Dl) 

da + db da + db 

1. Thermal field averaging 

Consider first the heat equation. Define an average temperature in the film, 6f{z), where 
we use the subscript F to indicate a quantity in the film suitably averaged over a period 
of the structure (averaging will be different for different quantities), under the assumption 
that the distance over which this averaged quantity varies significantly is much greater than 
the period da + db. Since the temperature field is continuous at the boundaries between the 
layers, Op{z) = 6a{z) = Ob{z). To obtain an averaged heat equation for the propagation of 
6f, begin with Eq. El here rewritten in a more convenient form, 

^u^C, e,{z) - ^ [k,^"^ = -luC, , (D2) 
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where q = a,b indicates a quantity evaluated in layer a or b, respectively. Averaging the 
first and last terms over a period of the structure is trivial. The second term requires more 
care. Noting that the continuity of the heat flux requires that kadOa/dz = kfydOh/dz = 
[k d9/dz)scvgi implicitly defining the averaged thermal conductivity kp by writing the average 
heat flux as 



\ / avg \ / avg 



and solving for kp, we find 



k-p" = (A;-i)avg . (D3) 
We can then write the averaged Eq. ID2I in a form analogous to Eq. 



iujOpiz) - kf-^-it = -iujPf , (D4) 



dz 



where the averaged film properties are 



Cf = (C) 



avg 



Kf = kp/CF 

where we used Eq. ^Jfor f3. With these averaged quantities in place of those of the uniform 
film, i.e. taking Xj — > Xp, we can immediately transcribe all the previous results for the 
temperature field in the uniform film without further analysis. It is also convenient to define 
a thermal diffusion time for the averaged film of thickness /, 

tf = P/kf . (D6) 



2. Elastic fields in a multilayer 

Averaging of the elastic properties is more straightforward. The zeroth-order elastic 
fields already taken as invariant through the region of interest, i.e. the in-plane dilation 
Eq = {eq^xx + ^o,to)/2 and the axial stress do = (Tq^zz, remain invariant in the multilayer, so 
they are obviously equal to their average. 

The calculation of the remaining components of the elastic field then follows exactly 
as given in appendix El so that the correct result for the fields in material g = a, 6 in a 
modulated film can be obtained from the corresponding expression for the zeroth-order field 
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in a uniform film j = / by replacing j — > q. For example, for the case Eq specified and 
o"o = we simply have, analogously to Eq. IA6I 

where 



A —A — 1 4 — ^'^g 



Eg 

Bo^xx,q = Bo^yy^q = , BQ^zz,q = . (D8) 

The same approach provides the results for ao,^,^, feo,M,g; ^o.m,.?; Sq, Sq, from the corrsponding 
expressions in Eqs. IA141 lAll IA81 and IA15[ respectively. 

For the power stored in the film, analogous to Eq. IA91 we must average the energy stored 
in the components of the multilayer, 

U. = (^) . (D9) 

By similar arguments as were applied to the zeroth-order fields, the first-order fields 
analogous to those obtained for a uniform film in Appendix [HI can be obtained from the 
corresponding expression for the first-order field in a uniform film J = / by replacing j — >■ q. 
For example, analogously to Eqs. \C6\ 

£l,ii,qiz) = Ai^ii^qaq9F{z), (Tl^ii^q{z) = B i ^qtt q9 p {z) , (DIO) 

where 

1 + z/, 



A, —A, — n 4, — ^ ' 

■'^l,xx,q -'^l,yy,q -'^l,zz,q ^ ^ 

Eq 

Bl^xx,q = Bi^yy q = —- , Bi :^z^q = . (DH) 

i Uq 

The combinations of these parameters used in calculating the dissipated power, Dg and dg, 
analogous to Eqs. IC8I and IC9t are obtained similarly. 



3. Averaging the dissipated power 

To find the averaged dissipated power, start with Eq. |2ZI Noting that the temperature 
is continuous and slowly varying over a period of the structure, we can write 

+ {DsSo + ds(Jo)as lm[-es{z)]dz^ . (D12) 
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Note that in making this approximation, we exclude cases where the thermal diffusion length 
approaches the layer period, but do allow the thermal length to be less than to the total 
thickness of the multilayer. This is not a very restrictive assumption: for typical mirror films 
of ~ 20 layer pairs, frequencies up to ~ 10^ above the dissipation peak are allowed (see Eq. 

. Comparing with Eq. |23 we see that any result for a uniform film can be transformed 
into the corresponding result for the averaged film by replacing 

{Df£o + dfao)af {{Deo + dao)a)^^^ , (D13) 

where the elastic quantities required are given in Eqs. IC8I and IC91 and the averaging oper- 
ation is defined in Eq. ID 11 we again replace the thermal properties with the averaged ones 
given in Eqs. ID3I and lD5l i.e. Xf — > Xp. Since the dissipated power is the key quantity 
from which all the end results of this paper devolve, only straightforward substitution and 
algebraic manipulation are required to obtain those results for the averaged film. The results 
so obtained for the dissipation factor and spectral density of thermoelastic noise Sx{f) are 
given in Eqs. EH and [721 respectively. 
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